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THERMOCHEMISTRY  AND  REACTIVITY  OF  METALS  ENGAGED  IN 
CHEMIIONIZATION 


Principal  Investigator:  Peter  B.  Armentrout  (armentrout@chem.utah.edu) 

Institution:  University  of  Utah,  Department  of  Chemistry 

Salt  Lake  City,  UT  841 12-0850 

Funding  Opportunity  Announcement  Number:  BAA-AFOSR-20 14-0001;  CFDA  12.800 

Air  Force  Office  of  Scientific  Research;  Energy,  Power  and  Propulsion  Sciences  (RTE) 
Molecular  Dynamics  and  Theoretical  Chemistry  (RTE-4) 

Program  Officer:  :  Dr.  Michael  R.  Berman  AFOSR/RTE  (michael.berman@us.af.mil) 

I.  Statement  of  Objectives 

The  thermochemistry  and  reactivity  of  oxidation  reactions  of  lanthanides  have 
implications  for  understanding  chemiionization  reactions  of  these  elemental  metals.  The  Air 
Force  has  explored  the  possibility  of  using  the  chemiionization  reactions  of  atomic  metal  ions 
with  oxygen  atoms  available  in  the  ionosphere  as  a  means  to  form  enhanced  plasma  densities  in 
the  atmosphere.  Such  man-made  plasmas  are  one  potential  means  to  mitigate  scintillation  effects 
that  interfere  with  satellite  communications.  Underscoring  the  need  for  fundamental  information 
regarding  such  reactions  are  recent  atmospheric  release  studies  involving  samarium  metal,  which 
failed  to  proceed  as  expected.  To  understand  this  failure,  it  became  clear  that  the  literature  data 
regarding  the  thermodynamics,  reactivity,  and  dynamics  of  samarium  oxidation  was  incomplete. 

To  acquire  thermodynamic  information  on  the  oxidation  of  samarium  and  other 
lanthanides,  reactions  of  the  atomic  metal  cations  and  their  oxides  with  atmospheric  gases  (O2, 
CO,  CO2,  NO,  H2O)  will  be  studied  as  a  function  of  kinetic  energy  in  a  guided  ion  beam  tandem 
mass  spectrometer  (GIBMS).  Such  reactions  are  of  direct  interest  in  understanding  the  chemistry 
that  might  occur  when  such  lanthanides  are  exposed  to  atmospheric  air  at  high  temperature. 
Analysis  of  such  data  will  provide  quantitative  bond  energies  of  a  variety  of  species  yielding  a 
better  fundamental  understanding  of  their  chemical  and  physical  properties.  Such 
thermodynamic  data  will  also  provide  better  predictability  regarding  the  reactivity  of  lanthanides, 
which  is  then  augmented  by  the  direct  kinetic  data  obtained  for  the  specific  reactions  examined. 
Both  the  kinetic  and  thermodynamic  data  will  be  useful  in  modeling  the  chemistry  of  these 
species. 

Accurate  quantum  chemical  calculations  on  heavy  metal  species  are  particularly 
problematic  because  of  the  need  to  implement  relativistic  quantum  chemistry  under  conditions  in 
which  electron  correlation  is  also  of  great  importance.  The  challenges  of  these  calculations  are 
exacerbated  by  a  shortage  of  accurate  information  on  such  heavy  metal  molecules,  which  makes 
validation  of  the  approximations  employed  difficult.  One  aim  of  the  present  work  was  to 
develop  a  base  of  accurate  thermochemical  data,  thereby  providing  benchmarks  against  which 
quantum  chemical  methods  may  be  tested. 

Accompanying  its  scientific  merit  is  the  fact  that  this  work  involved  the  education  of 
graduate  and  undergraduate  students  in  lanthanide  science  and  experiments. 
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II.  Research  Effort 
A.  Background 

Satellite  communications  (e.g.,  global  positioning  systems,  GPS,  and  telecommuncations) 
are  of  clear  importance  in  both  the  civilian  and  military  sectors,  however,  scintillation  effects  can 
interfere  with  radio  wave  propagation  leading  to  intermittent  or  complete  loss  of  communication. 
Scintillation  can  be  caused  by  natural  irregularities  in  the  ionosphere.  In  critical  applications,  it 
may  be  desirable  to  mitigate  these  scintillation  effects  by  directly  inducing  man-made  enhanced 
plasma  densities  to  slow  the  growth  of  the  ionospheric  instabilities.  The  Air  Force  has  explored 
the  possibility  of  creating  such  artificial  plasma  clouds  by  increasing  the  electron  density  using 
chemical  release.  Materials  such  as  barium,  strontium,  xenon,  lithium,  and  cesium  have  been 
studied,  but  require  solar  photons  or  particle  collisions  for  ionization,  hence  limiting  their 
efficiency.  Other  processes  that  hold  higher  promise  in  this  regard  are  chemiionization  processes, 
and  in  particular,  reaction  1  where  M  is  an  elemental  metal. 

M  +  O  ->  MO+  +  e~  AH  <  0  kJ/mol  (1) 

If  reaction  1  is  exothermic,  this  process  can  occur  spontaneously  by  utilizing  the  atomic 
oxygen  present  at  the  altitudes  of  interest.  Further,  the  undesirable  reverse  reaction  is  inefficient 
because  it  is  endothermic.  (This  simplistic  conclusion  is  mediated  by  the  fact  that  the  electron 
temperature  in  the  ionosphere  is  -1000  K,  such  that  some  appreciable  fraction  of  energetic 
electrons  can  undergo  the  reverse  reaction  if  reaction  1  is  not  exothermic  by  more  than  -0.4  eV.) 
Recently,  Air  Force  Research  Laboratory  (AFRL)  Metal  Oxide  Space  Cloud  (MOSC) 
experiments  have  been  conducted  to  test  the  efficiency  of  forming  an  artificial  plasma  cloud  by 
this  process  using  samarium,  M  =  Sm.  Samarium  was  chosen  for  these  initial  experiments 
because  the  literature  thermochemistry  indicated  the  reaction  was  exothermic  by  about  0.3  eV,1 
and  samarium  has  a  relatively  low  boiling  point  (2067  K),  which  would  allow  efficient 
vaporization  of  this  material  using  a  titanium-boron  thermite  mixture  that  reaches  temperatures 
of  about  3500  K.  Predictions  of  the  test  suggested  that  only  a  few  kg  of  metal  dispersed  in  the 
ionosphere  near  an  altitude  of  170  km  could  create  an  enhanced  plasma  cloud  over  distances 
extending  as  far  as  100  km.  Unfortunately,  the  plasma  density  produced  was  about  two  orders  of 
magnitude  smaller  than  expected,  and  the  cloud  shape  and  color  all  differed  significantly  from 
pre-launch  predictions. 

To  understand  these  phenomena,  the  plasma  chemistry  research  group  of  Dr.  Albert 
Viggiano  at  AFRL  was  contacted  to  explore  the  chemistry  of  samarium  more  thoroughly.  Using 
their  selected  ion  flow  tube  (SIFT)  apparatus  equipped  with  an 
electrospray  ionization  (ESI)  source,  they  determined  the  rates  and 
temperature  dependences  of  several  reactions  of  Sm+  and  SmO+ 
near  room  temperature,  finding  that  N20,  N02,  02,  and  S02  all 
react  with  Sm+  to  form  SmO+.2  Examination  of  the 
thermochemistry  in  Table  1  shows  that  these  observations  suggest 
that  D(Sm+-0)  >  5.66  eV  =  D(OS-O).  In  contrast,  the  SIFT  studies 
also  found  that  C02  and  NO  were  unreactive  with  Sm+.  Because 
D(OC-O)  is  lower  than  5.66  eV  (Table  1),  the  failure  to  react  with 
C02  at  room  temperature  cannot  be  a  consequence  of  unfavorable 
thermodynamics,  but  rather  must  result  from  kinetic  effects.  In 
addition,  the  SIFT  studies  determined  that  SmO+  reacted  with  N20, 

N02,  02,  C02,  S02,  and  NO  exclusively  by  three-body  association. 


Table  1 


Oxide  Bond  Energies 
Species  Dq(M-O),  eV 


n2o 

1.672  +  0.004 

no2 

3.116  +  0.009 

h2o 

5.03410.001 

02 

5.115  ±0.001 

co2 

5.453  10.002 

so2 

5.6610.02 

NO 

6.507  1 0.002 

CO 

11.108  ±0.005 

3 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Unfortunately,  these  data  yielded  no  specific  insight  into  the  MOSC  results,  hence,  the  Pis 
laboratory  was  contacted  to  see  if  we  could  elucidate  the  thermochemistry  of  SmO  and  SmO+ 
more  thoroughly. 

As  described  further  below,  the  guided  ion  beam  tandem  mass  spectrometry  (GIBMS) 
approach  developed  in  the  Pis  laboratory  has  proven  to  be  a  versatile  experimental  method  of 
determining  quantitative  thermodynamic  information.  In  this  particular  case,  the  PI  and  his 
graduate  student,  Richard  M  Cox,  studied  several  reactions  that  provided  specific  information 
regarding  the  bond  energy  of  SmO+,  with  the  results  having  recently  been  published.2  These  data 
indicate  a  substantially  lower  bond  energy  than  suggested  in  the  literature,  sufficiently  so  that 
reaction  1  now  appears  to  be  exothermic  by  only  0.08  ±  0.07  eV.  Simulations  of  the  MOSC  test 
using  this  thermodynamic  information  are  in  agreement  with  the  experimental  results  from  this 
test,  thereby  confirming  the  utility  of  the  fundamental  thermodynamic  data  provided  by  the  Pis 
laboratory.  Continued  work  on  these  systems  has  extended  our  results  to  that  of  gadolinium 
oxide  (Gd),  as  studied  by  graduate  student,  Maria  Demireva. 

B.  Experimental  methods 

Overview.  The  experimental  technique  we  use  is  GIBMS.3  This  device  uses  mass 
spectrometry  to  select  an  ionic  reactant,  which  reacts  under  single  collision  conditions  with 
neutral  molecules  in  a  radiofrequency  (rf)  octopole  ion  beam  guide,  and  then  analyzes  the 
products  using  mass  spectrometry,  with  single  molecule  sensitivity.  The  Utah  GIBMS  instrument 
allows  the  study  of  the  kinetic  energy  dependence  of  ion-molecule  reactions  over  a  wide  (four 
order  of  magnitude)  energy  range.  This  enables  the  measurement  of  the  thermodynamics  of 
metal-ligand  interactions  using  two  types  of  processes:  exchange  reactions  2, 

M+  +  LR  ->  ML+  +  R  (2) 

and  collision-induced  dissociation  (CID),  reactions  3. 

MLX+  +  Xe  ->  MLx_i+  +  L  +  Xe  (3) 

In  the  first  case,  if  the  reaction  is  endothermic,  the  threshold  for  reaction,  Eq,  can  be  measured 
and  related  to  the  bond  dissociation  energy  (BDE)  of  M+-L  by  Do(M+-L)  =  Do(LR)  -  Eq,  where 
Do(LR)  is  taken  from  the  literature.  In  the  CID  process,  which  is  intrinsically  endothermic,  Eq 
equals  D0(Lx_iM+-L)  as  long  as  there  are  no  barriers  to  the  reverse  process,  which  is  generally  the 
case  because  of  the  long-range  attractive  character  of  ions  with  polarizable  and  polar  molecules.4 
Accurate  determination  of  metal-ligand  BDEs  using  this  approach  requires  that  several  factors  be 
taken  into  account,  as  detailed  below.5  6 

Ion  generation.  The  ion  source  is  a  critical  feature  for  the  quantitative  study  of  the 
energetics  of  ion-molecule  reactions.  Unless  the  internal  energy  of  the  ions  is  adequately 
controlled,  our  knowledge  of  the  energy  available  to  reaction  and  thus  the  precision  and  accuracy 
of  any  thermodynamic  measurement  is  compromised.  For  the  studies  proposed  here,  this  control 
is  achieved  using  a  dc  discharge/flow  tube  (DC/FT)  source.7  13  This  source  has  been  used  to 
generate  atomic  metal  and  metal  oxide  ions  for  studies  of  reactions  2  and  3  with  most  of  the 
transition  metals  (TMs),  measuring  a  wide  variety  of  covalent  metal-ligand  bonds.14"15  It  is  also 
the  source  used  in  the  samarium  and  preliminary  gadolinium  chemistry  outlined  below,  as  well 
as  recent  work  on  the  actinide  thorium.16  In  this  source,  argon  (about  10%  in  He)  is  ionized  by  a 
kV  dc  potential  and  then  collides  with  the  cathode  composed  of  or  containing  the  material  of 
interest,  essentially  sputtering  the  desired  metal  cations.  Ions  generated  in  this  glow  discharge 
can  then  undergo  further  reaction  (such  as  oxidation  or  condensation)  in  the  meter-long  flow 
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tube.  In  addition,  the  ions  are  thermalized  by  undergoing  about  105  collisions  with  the  inert  flow 
gases.  For  atomic  transition  metal  cations,  we  have  demonstrated  that  the  distribution  of 
electronic  states  produced  by  the  DC/FT  source  is  thermal  or  nearly  so.  For  instance,  by 
comparing  results  with  those  from  a  surface  ionization  source  (where  the  source  temperature  is 
known),  the  DC/FT  source  has  been  found  to  generate  Sc+,17  Fe+,18  Co+,19  Ni+,20  Ru+,21  Rh+,21 
and  Pd+  21  ions  with  an  average  electronic  temperature  of  700  ±  400  K,  and  Y+,  Zr+,  Nb+,  and 
Mo  ions  with  an  average  electronic  temperature  of  300  ±100  K.  Therefore,  we  conservatively 
quote  the  DC/FT  source  temperature  as  700  ±  400  K.  For  lanthanide  cations,  this  means  that  the 
dominant  species  occupied  is  the  ground  state  term. 

As  detailed  below,  the  DC/FT  source  can  produce  species  such  as  M+,  MO+,  and  MC>2+ 
for  M  =  Sm  and  Gd.  We  have  also  demonstrated  that  termolecular  reactions  can  attach  one  or 
many  ligands  of  all  sorts  (Ar,  CO,  N2,  H2O,  CO2,  amino  acids,  through  crown  ethers)  to  bare 
metal  and  oxygenated  ions  throughout  the  periodic  table,  '  ’  ’  ’  thereby  generating  a  wide 
range  of  MLX+  species,  e.g.,  M+(02),  M+(C02)  and  OM+(CO).  Using  reactions  3,  quantitative 
measurements  of  noncovalent  metal-ligand  interactions  have  been  examined  for  alkalis,  alkaline 
earth,  and  transition  metal  ions  with  a  wide  range  of  ligands.14  15  25  Below,  such  studies  are 
illustrated  by  determination  of  the  bond  energies  of  Sm+-C02,  OSm+-CO,  and  Gd+(02). 

Guided  ion  beam  mass  spectrometry.  Once  the  ions  leave  the  source,  they  pass  through 
a  differentially  pumped  chamber  where  they  are  focused  into  a  beam  and  mass  analyzed.  They 
then  enter  another  chamber,  where  they  are  decelerated  to  a  desired  kinetic  energy  and  injected 
into  an  rf  octopole  beam  guide,26  which  acts  as  a  radial  ion  trap.  The  octopole  passes  through  a 
reaction  cell  containing  a  neutral  reactant  gas  maintained  at  sufficiently  low  pressure  that  single¬ 
collision  conditions  prevail.  The  pressure  dependence  of  all  cross  sections  is  measured  to  ensure 
that  rigorous  single-collision  conditions  are  used  in  the  data  analysis.  Product  and  remaining 
reactant  ions  drift  to  the  end  of  the  octopole,  where  they  are  focused  into  a  final  differentially 
pumped  chamber.  (In  one  instrument,  the  octopole  region  contains  two  octopoles,  which  allows 
additional  dynamics  experiments  to  be  conducted.  )  Here,  the  ions  are  separated  with  a 
quadrupole  mass  filter  and  detected  using  a  secondary  electron  scintillation  (Daly)  detector”  and 
standard  ion  pulse  counting  techniques  that  provide  essentially  unit  detection  capability  up  to 
high  masses.  Data  collection  is  under  computer  control,  which  allows  extensive  signal  averaging. 
At  each  relative  kinetic  energy,  E ,  the  intensities  of  transmitted  reactant  and  product  ions  are 
converted  to  absolute  reaction  cross  sections  for  each  product  channel,  a(E).  These  cross 
sections  have  absolute  uncertainties  of  ±20%.  Relative  cross  sections  are  accurate  to  about  ±5%. 

There  are  several  advantages  to  using  guided  ion  beam  techniques  to  study  ion-molecule 
reactions.  First,  the  range  of  energies  available  to  the  instrument  is  extensive,  four  orders  of 
magnitude  ranging  from  thermal  to  a  thousand  electron  volts.  This  is  the  primary  distinction 
between  this  method  of  examining  gas-phase  ion-molecule  chemistry  and  other  instrumentation, 
which  is  generally  limited  to  thermal  energies.  Second,  the  absolute  kinetic  energy  scale  can  be 
determined  accurately  using  the  octopole  as  a  highly  efficient  retarding  energy  analyzer.  Third, 
the  octopole  ion  beam  guide  allows  for  highly  efficient  product  collection.  ’  ’  Fourth,  because 
the  octopole  extends  well  beyond  the  ends  of  the  gas  cell  (unlike  triple  quadrupole  collision 
regions),  the  collision  energy  is  well-controlled  and  the  1/r6  effective  potential  is  much  less 
perturbing  than  the  1/r2  potential  in  quadrupoles. 

Data  analysis.  The  ability  to  examine  reactions  over  a  very  wide  range  of  kinetic 
energies  is  the  distinguishing  and  truly  powerful  ability  of  guided  ion  beam  methods.  This 
capability  allows  endothermic  reactions  to  be  examined  routinely.  By  measuring  the  threshold 


5 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


for  reaction,  we  can  obtain  thermodynamic  information  on  the  reaction  products  for  a  wide 
variety  of  species.  This  generally  assumes  that  activation  barriers  in  excess  of  the  endothermicity 
are  not  present.  This  assumption  is  generally  valid  for  ion-molecule  reactions,  largely  because  of 
the  long-range  attractive  form  of  ion-molecule  interactions,  as  long  as  there  are  no  constraints 
resulting  from  spin  or  orbital  angular  momentum  conservation.5  (Interestingly,  the  samarium 
chemistry  provides  a  couple  of  examples  where  such  unusual  constraints  are  evident,  although 
not  yet  completely  understood,  as  discussed  further  below.)  For  CID  studies  involving 
heterolytic  bond  cleavages,  which  encompasses  many  of  the  systems  of  interest  here,  we  have 
shown  that  there  should  be  no  reverse  activation  barriers  because  of  electronic  considerations.4 

Our  methods  of  analysis  have  been  amply  described  in  the  literature,5"6  and  have  proven 
to  be  robust  for  a  large  number  of  species.  Briefly,  we  use  Eq.  4  to  reproduce  the  energy 
dependence  of  reaction  cross  sections. 

=  J(1  -e-^^’KE-ey-'de  (4) 

Here,  oo  is  an  energy  independent  scaling  factor,  E  is  the  relative  translational  energy  of  the  ion, 
Eo  is  the  reaction  threshold  at  0  K,  and  n  is  an  adjustable  fitting  parameter  that  describes  the 
efficiency  of  the  energy  transfer  upon  collision.”  The  summation  is  over  the  rovibrational  and 
electronic  states  of  the  reactants  having  excitation  energies,  Eh  and  populations,  g,-,  where  Eg,  = 
1 .  e  represents  the  energy  deposited  into  the  ion  upon  collision  with  the  neutral  reagent,  and  r  is 
the  average  timescale  of  the  experiment  (ion  time-of-flight  from  the  collision  cell  to  the 
quadrupole  mass  spectrometer,  either  ~0.1  or  0.5  ms  for  the  single  and  dual  octopole 
instruments,  respectively).  The  rate  constant,  k{£  +  £))  =  k(E*),  is  the  Rice-Ramsperger-Kassel- 
Marcus  (RRKM)  unimolecular  dissociation  rate  constant31’32  as  calculated  using  Eq.  5, 

k(E*)  =  sNfvr  (E*  -  Eo)  /  hfrr  (E*)  (5) 

where  s  is  the  reaction  degeneracy  calculated  from  the  ratio  of  rotational  symmetry  numbers  of 
the  reactants  with  respect  to  the  products,  Nfvr(E*  -  E0)  is  the  sum  of  the  rovibrational  states  of 
the  transition  state  (TS)  at  an  energy  E*  -  Eo  above  the  threshold,  Eo,  and  Pvr(E*)  is  the  density 
of  rovibrational  states  for  the  energized  molecule  (EM)  at  the  energy  available,  E*.  When  the 
rate  constant  is  much  faster  than  the  average  experimental  timescale  (which  is  generally  true  for 
small  molecular  systems),  Eq.  4  reduces  to  Eq.  6, 

o(E)=obZgi(E  +  Ei-Eo)n/E  (6) 

a  simple  modified  line-of-centers  form.  This  form  is  often  sufficient  to  analyze  the  cross  sections 
for  exchange  reactions  2.  After  convolution  over  the  kinetic  energy  distributions  of  both 
reactants,  Eqs.  4  and  6  have  been  shown  to  accurately  describe  the  kinetic  energy  dependence  of 
many  bimolecular  reactions  and  CID  experiments  (see  examples  below  in  Figures  2  and  3). 

This  model  for  the  kinetic  energy  dependence  of  ion-molecule  reactions  includes  explicit 

o  oi  oc 

consideration  of  the  internal  energy  of  both  reactants,  '  ’  and  has  been  shown  to  be  directly 
compatible  with  direct  measurements  of  the  energy  transfer  function.  For  CID  of  larger 
systems,  the  model  includes  the  possibility  that  the  collisionally  excited  complex  does  not 
dissociate  during  the  experimental  flight  time,  which  can  lead  to  a  delayed  onset  for  the  CID 
threshold,  a  kinetic  shift.  Our  model  includes  an  estimate  of  this  effect  by  incorporating 
statistical  rate  theory  for  the  unimolecular  rate  of  dissociation  of  an  energized  molecule.36"37 
Several  reasonable  treatments  of  the  assumed  energy  and  angular  momentum  distributions  of  the 
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collisionally  activated  molecules  were  considered.  This  statistical  approach  also  allows  a 

oo 

straightforward  introduction  of  competing  dissociation  channels,  for  which  results  demonstrate 
that  a  statistical  model  for  the  rotational  energy  distributions  gives  the  most  accurate 
thermodynamic  data.  Additional  work  shows  that  the  best  results  may  require  explicit  handling 

on 

of  internal  rotations.  These  statistical  analysis  methods  have  been  extensively  tested  and  further 
developed.40  45  These  studies  demonstrate  that  a  phase  space  theory  (PST)  approach  in  which 
rotational  and  orbital  angular  momentum  are  explicitly  conserved  and  allowed  to  interchange 
provides  the  best  reproduction  of  the  data,  however,  they  also  demonstrate  that  the 
thermodynamic  information  obtained  is  insensitive  to  the  assumptions  regarding  angular 
momentum.  Overall,  the  model  has  been  shown  to  reproduce  the  threshold  regions  of  a  variety  of 
endothermic  ion-molecule  reactions,  including  CID,  with  good  accuracy  in  the  resulting 
thermochemistry.  Recent  efforts  have  now  extended  our  “toolbox”  of  thermodynamic 

10 

methodology  from  simple  bimolecular  reactions  and  CID,  to  variants  involving  competitive 
and  sequential  dissociation  pathways,46  as  well  as  to  association44  and  ligand  exchange 
reactions.45  These  methods  therefore  allow  multiple  pathways  for  obtaining  and  verifying 
accurate  bond  energies,  as  well  as  free  energies  and  entropies.  Although  the  precision  of  such 
measurements  does  not  approach  spectroscopic  analogues,  it  is  nevertheless  still  quite  good, 
often  less  than  0.1  eV  and  occasionally  as  low  as  0.02  eV.47 

C.  Results 

Reactions  involving  samarium.  The  literature  regarding  the  BDE  of  SmO,  Do(SmO), 
has  been  reviewed  thoroughly  in  our  recent  paper.2  Data  from  several  high  temperature  studies 
yield  BDEs  with  a  range  of  5.5  ±  0.2  -  6.10  ±  0.03  eV48  50  and  have  been  analyzed  in  several 
compilations.5155  The  most  inclusive  and  conservative  value  is  5.88  ±  0.17  eV.52  The  ionization 
energy  of  SmO,  IE(SmO),  has  been  measured  using  electron  ionization  and  high  temperature 
mass  spectrometry  as  5.55  ±0.1  eV.56  These  two  pieces  of  information  alone  allow  the  enthalpy 
of  reaction  1,  AH0(1)  =  IE(MO)  -  D0(MO),  to  be  determined  as  5.55  -  5.88  =  -0.33  ±  0.20  eV. 
Combined  with  IE(Sm)  =  5.6437  eV,57  these  data  also  indicate  that  the  BDE  of  SmO+,  Do(Sm+- 
O),  equals  5.97  ±  0.20  eV  because  IE(MO)  -  D0(MO)  =  IE(M)  -  D0(MO+).  This  BDE  is 
consistent  with  the  observations  of  the  AFRL  laboratory  using  their  SIFT  apparatus  that  reaction 
of  Sm+  with  NoO,  N02,  (A,  and  S02  were  all  exothermic  to  form  SmO+,  as  well  as  the  lack  of 
reaction  with  NO  as  it  is  endothermic  (Table  1).“  This  BDE  also  indicates  that  reaction  with  C02 
is  exothermic,  whereas  no  reactivity  was  observed  at  thermal  energies  between  Sm+  and  C02  in 
the  AFRL  SIFT  studies.2 

In  our  GIBMS  studies  of  the  reactions  of  Sm+,  the  ions  were  generated  using  the  DC/FT 
source.  For  Sm+,  an  electronic  temperature  of  700  ±  400  K  means  that  21  -  68  %  of  the  ions  are 
in  their  ground  8Fj/2  level  with  27  -  34%  in  the  xFy2  level  at  0.04  eV,  and  4  -  21%  in  the  ’sF2/2 

co 

level  at  0.104  eV.  Overall,  the  average  electronic  energy  is  0.06  ±  0.05  eV,  an  uncertainty  that 
is  included  in  our  determination  of  the  threshold  energies. 

Our  results  for  reaction  of  Sm+  with  02,  S02,  C02,  and  NO  agree  with  the  AFRL  findings 
and  are  shown  in  Figure  1.  Figures  la  and  lb  show  that  the  reactions  of  Sm+  with  02  and  S02 
efficiently  form  SmO+  with  cross  sections  exhibiting  no  barrier  as  low  in  energy  as  we  can  go, 
suggesting  that  Do(SmO+)  >  D0(O2)  =  5.115  eV  and  Do(SmO+)  >  Do(OS-O)  =  5.66  eV.  In  the  02 
system,  the  cross  section  matches  the  predicted  Langevin-Gioumousis-Stevenson  (LGS)  collision 

in 

cross  section  for  ion-molecule  reactions,  suggesting  a  reaction  efficiency  of  100  ±  20%.  This 
efficiency  differs  from  the  AFRL  SIFT  result  of  about  50%,  which  agrees  with  a  previous 
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Fig.  1.  Cross  sections  for  reaction  of  Sm+  with  O2,  SO2,  CO2,  and  NO  as  a  function  of  kinetic 
energy  in  the  center-of-mass  (lower  axis)  and  laboratory  (upper  axis)  frames.  Arrows  indicate  the 
bond  energies  of  the  neutral  reactants.  Oscillations  in  the  cross  section  for  SmO+  in  parts  a  and  d 
are  artifacts  of  the  ion  focusing.  The  Sm02+  product  in  part  a  shows  a  linear  dependence  on  the 
pressure  of  O2. 


literature  report  (also  SIFT  work),59  a  discrepancy  that  is  not  understood  although  it  could  be 
related  to  the  different  ion  sources  used.  For  the  SO2  reaction,  the  SmO+  cross  section  parallels 
the  collision  cross  section  (calculated  here  using  the  Su-Chesnavich  semi-classical  trajectory 
method  to  account  for  the  dipole  of  SO2)60  at  low  energies  and  has  an  efficiency  of  22  ±  4%, 
similar  to  the  SIFT  results,  29  ±  9%. 

In  the  O2  system,  the  Sm02+  product  is  observed  at  higher  energies,  but  this  cross  section 
depends  linearly  on  the  O2  pressure  indicating  that  it  is  formed  by  reaction  of  the  SmO+  product 
with  another  molecule  of  O2.  In  the  SO2  system,  several  minor  products  are  observed.  Figures  lc 
and  Id  clearly  show  that  the  reactions  of  Sm+  with  CO2  and  NO  exhibit  barriers  (a  similar  result 
is  also  obtained  for  oxidation  of  Sm+  by  OCS),  consistent  with  the  failure  to  observe  any  reaction 
at  thermal  energies  in  the  SIFT  studies.  Our  analysis  of  these  cross  sections  using  the  methods 
discussed  above  yield  thresholds  for  these  reactions  of  1.84  ±  0.12  eV  and  2.59  ±  0.16  eV, 
respectively.  If  these  thresholds  corresponded  to  the  asymptotic  energies  of  the  products,  they 
could  be  combined  with  the  CO2  and  NO  bond  energies,  Table  1,  to  yield  D0(Sm+-O)  =  3.61  ± 
0.12  and  3.92  ±  0.17  eV,  respectively.  These  values  cannot  correspond  to  the  true  BDE  otherwise 
the  oxidation  of  Sm+  by  O2  and  SO2  would  be  endothermic.  Thus,  these  thresholds  must  actually 
correspond  to  barriers  along  the  potential  energy  surfaces,  a  result  that  can  be  explicitly 
demonstrated  by  examining  the  reverse  reaction,  as  discussed  further  below.  These  barriers  are 
actually  fairly  unusual  in  ion-molecule  chemistry  and  continue  to  be  explored.  Their  elucidation 
should  reveal  interesting  electronic  properties  of  lanthanide  chemistry. 


8 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


We  also  examined  several  reactions  not  previously  studied  by  the  AFRL  group,  including 
Sm+  +  CO,  SmO+  +  Xe,  and  SmO+  +  CF.  Our  results  for  these  processes  are  shown  in  Figure  2, 
where  the  modeling  of  the  data  is  also  shown.  (Reaction  of  SmO+  +  O2  also  forms  Sm02+  in  an 
endothermic  process  that  allows  D0(OSm+-O)  =  1.14  ±  0.15  eV  to  be  measured.2)  For  reaction 
with  CO,  the  threshold  for  SmO+  production  is  5.49  ±  0.12  eV,  which  combined  with  Do(CO)  in 
Table  1,  yields  Do(SmO+)  =  5.62  ±  0.15  eV,  which  is  consistent  with  the  exothermic  reaction  of 
Sm+  +  O2  and  SO2  within  experimental  uncertainty.  In  the  collision-induced  dissociation  (CID) 
reactions  with  Xe  and  O2,  the  thresholds  measured  are  5.67  +  0.16  and  5.78  ±  0.09  eV, 
respectively.  These  values  should  directly  correspond  to  Do(SmO+),  although  conservatively  they 
are  upper  limits  as  the  translational  to  internal  energy  transfer  can  be  inefficient  for  strongly 
bound  species  like  SmO+.  Notably,  these  three  BDEs  are  within  experimental  uncertainty  of  each 
other  and  have  a  weighted  average  of  D0(Sm+-O)  =  5.72s  +  0.07  eV  (one  standard  deviation).2 
This  value  also  agrees  with  the  best  literature  value,  5.97  ±  0.20  eV,  within  combined 
experimental  uncertainties,  but  is  more  precise  and  lower  by  0.25  eV. 
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Fig.  2.  Cross  sections  for  reaction  of  Sm+  with  CO  (part  a)  and  SmO+  +  Xe  and  O2  (part  b)  as  a 
function  of  kinetic  energy  in  the  center-of-mass  (lower  x-axis)  and  laboratory  (upper  x-axis) 
frames  (for  Xe  in  part  b).  Dashed  lines  are  models  of  the  cross  sections,  and  solid  lines  are  the 
same  models  convoluted  over  the  kinetic  and  internal  energy  distributions  of  the  reactants.  In 
part  a,  the  arrow  indicates  the  bond  energy  of  the  neutral  reactant. 
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An  important  test  of  this  value  was  also  reported  in  our  recent  work.”  The  group  of  M. 
Heaven  (Emory  U.)  measured  IE(SmO)  using  resonantly  enhanced  two-photon  ionization 
measurements  (REMPI)  and  pulsed-field  ionization  zero  kinetic  energy  (PFI-ZEKE) 
photoelectron  spectroscopy.  Their  final  result,  46318  ±  5  cm"1  (5.7427  ±  0.0006  eV),  is  0.19  eV 
higher  than  the  literature  value  of  5.55  ±0.1  eV.56  In  previous  work,61  63  the  IEs  reported  by 
Rauh  and  Ackermann56'64  for  ZrO,  HfO,  and  TaO  have  also  been  found  to  be  too  low  by  even 
larger  amounts.  This  systematic  difference  may  occur  because  the  population  of  excited  states  at 
the  elevated  temperatures  used  in  these  early  experiments  was  not  adequately  accounted  for. 
Combined  with  the  literature  value  for  D0(SmO)  =  5.88  ±  0.17  eV  and  IE(Sm)  =  5.6437  ±  0.0006 
eV,  this  refined  IE(SmO)  yields  Do(SmO+)  =  5.78  ±  0.17  eV,  in  agreement  with  the  GIBMS 
value.  Conversely,  combining  these  precise  IEs  with  our  value  for  D0(SmO+)  yields  a  refined 
neutral  BDE  of  D0(SmO)  =  5.83  ±  0.07  eV. 

Finally,  we  return  to  the  enthalpy  for  reaction  1.  According  to  our  measurements,  the 
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SmO+  bond  energy  is  5.72s  ±  0.07  eV  (weighted  average  of  three  values).  Given  this  value  and 
IE(Sm),  the  chemiionization  reaction  1  is  exothermic  by  0.08  ±  0.07  eV,  0.25  eV  lower  than 
previously  thought.  Recent  models  of  the  MOSC  test  using  this  thermochemistry  have  nicely 
accounted  for  the  observations  made  in  the  atmospheric  tests.  This  result  is  also  consistent  with 
the  original  observations  of  both  Fite  and  coworkers65  and  Dyke  and  coworkers66  that  the 
chemiionization  process  with  samarium  is  observed  with  relatively  small  probabilities. 

As  noted  above,  the  observation  of  an 
appreciable  barrier  for  the  reactions  with  CO2, 
OCS,  and  NO  are  interesting  results  that  need 
to  be  explored  more  thoroughly,  both 
experimentally  and  computationally.  In  the  case 
of  the  CO2  reaction,  we  have  also  examined  the 
reverse  reaction,  SmO+  +  CO  -»  Sm+  +  CO2,  as 
shown  in  Figure  3.  Comparison  of  this  cross 
section  with  that  for  CID  with  Xe  or  O2  (Figure 
2b)  shows  that  the  sharp  rise  in  the  formation  of 
Sm+  at  energies  above  ~6  eV  can  be  attributed 
to  CID.  Therefore,  the  reactivity  observed  at 
low  energies  must  correspond  to  formation  of 
the  CO2  neutral  product.  This  result 
unambiguously  shows  that  there  is  a  barrier  in 
excess  of  the  product  asymptotes  for  both  the 
forward  and  reverse  reactions.  The  onset  of  this 
reverse  reaction  has  a  threshold  energy 
consistent  with  the  measured  barrier  height  of 
the  forward  reaction,  1.84  ±  0.12  eV,  combined 
with  the  calculated  exothermicity  of  this 
reaction,  D0(Sm+-O)  -  D0(OC-O)  =  5.725 
-  5.453  eV  =  0.27  ±  0.07  eV,  i.e.,  a  reverse 
barrier  height  of  2.1 1  ±0. 14  eV,  Figure  3. 

To  further  explore  this  potential 
energy  surface,  we  have  also  examined 
the  CID  reactions  of  OSm+(CO)  and 
Sm+(C02),  species  that  are  easily  formed 
in  the  flow  tube  source  by  attaching  the 
CO  or  CO2  molecules  to  SmO+  and  Sm+ 
using  three-body  association  reactions. 
These  experiments  establish  well  depths 
of  D0(OSm+-CO)  =  0.99  ±  0.07  eV  and 
D0(Sm+-CO2)  =  0.43  ±  0.03  eV.  Thus, 
five  points  along  the  Sm+  +  CO2  potential 
energy  surface  (reactants,  products,  two 
intermediates,  and  the  barrier  height)  are 
experimentally  established.  These  are 
shown  in  Figure  4  along  with  similar 
information  on  other  oxidation  reactions 


Fig.  3.  Cross  section  for  reaction  of  SmO+ 
with  CO  as  a  function  of  kinetic  energy  in 
the  center-of-mass  (lower  x-axis)  and 
laboratory  (upper  x-axis)  frames.  Dashed 
lines  are  models  of  the  cross  sections,  and 
the  solid  line  is  the  sum  of  these  models 
convoluted  over  the  kinetic  and  internal 
energy  distributions  of  the  reactants. 


Fig.  4.  Potential  energy  surfaces  for  oxidation  of 
Sm+  with  O2,  CO2,  SO2,  NO,  and  CO.  Small 
numbers  indicate  experimentally  measured  relative 
energies  in  eV.  The  location  of  the  excited  13/2 
state  of  Sm+  is  shown  by  the  dashed  line. 


10 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Energy  (eV,  Lab) 
10 


15 


studied  here.  Theoretical  exploration  of  these  potential  energy  surfaces,  including  consideration 
of  spin-orbit  effects,  is  now  ongoing  to  complete  our  examination  of  these  interesting  systems. 
Our  hypothesis  for  the  origins  of  the  barrier  are  also  indicated  in  Figure  4.  Namely,  formation  of 
a  strong  SmO+  bond  requires  at  least  two  electrons  in  the  valence  shell  (non  4f).67  The  lowest 
energy  state  of  Sm+  that  meets  this  requirement  is  a  13/2  state  having  a  4f55d16s1  configuration, 
3.29  eV  above  the  8Fi/2(4f66s1)  ground  state.  The  barriers  are  conjectured  to  be  a  result  of  the 
diabatic  curve  crossings  between  the  surfaces  evolving  from  these  two  states.  An  interesting 
question  remains  why  the  surfaces  for  the  02  and  S02  reactants  have  small  barriers  (lower  than 
the  reactants)  whereas  C02  and  NO  have  large  barriers,  which  may  be  related  to  the  spin  states  of 
the  neutral  products,  see  Figure  4. 

Another  example  of  this  type  of  study  is  shown  in  Figure  5.  Here  we  show  preliminary 
data  for  the  gadolinium  system  in  which  Gd02+  is  formed  in  the  flow  tube  by  addition  of  02, 
such  that  conceivably  either  Gd02+  (the  dioxide)  or  Gd+(02)  (the  adduct)  could  be  formed.  The 
CID  cross  section  with  Xe  clearly  exhibits  the 
formation  of  both  Gd+  +  02  and  GdO+  +  O,  with 
the  latter  channel  having  both  a  low  and  high 
energy  feature.  (Although  the  low  energy  feature  is 
noisy,  its  magnitude  is  decidedly  above  the 
baseline  noise  at  the  lowest  energies  below  0.3 
eV.)  What  is  not  shown  is  that  the  magnitude  (but 
not  the  kinetic  energy  dependence)  of  the  Gd+  and 
low-energy  GdO+  cross  sections  vary  appreciably 
in  magnitude  (by  over  a  factor  of  10)  as  the  source 
conditions  are  changed,  whereas  the  higher  energy 
GdO+  feature  remains  static.  This  behavior 
combined  with  the  relative  thresholds  indicates 
that  both  low  energy  features  arise  from 
dissociation  of  a  minor  amount  of  the  Gd+(02) 
adduct,  whereas  the  high  energy  GdO+  feature 
comes  from  the  Gd02+  dioxide,  which  dominates 
the  population  of  these  reactants.  Note  that  the 

low-energy  GdO+  cross  section  is  generally  smaller  than  the  Gd+  cross  section  but  does  appear  at 
slightly  lower  energies.  This  behavior  is  consistent  with  facile  loss  of  02  from  the  Gd+(02) 
adduct  whereas  loss  of  O  requires  surmounting  the  insertion  barrier  leading  to  the  Gd02+  dioxide 
and  therefore  is  entropically  disfavored.  Analyses  of  these  data  along  with  the  reactions  Gd+  +  02 
— »  GdO+  +  O  (which  like  Sm+  in  Figure  1  occurs  efficiently  with  no  barrier)  and  GdO+  +  02 
Gd02+  +  O  (which  is  endothermic)  means  that  we  can  measure  the  depths  of  the  potential  wells 
for  both  the  adduct  and  the  dioxide  as  well  as  the  barrier  height  separating  them  (from  the 
threshold  for  the  low-energy  GdO+  feature  in  Figure  5).  This  latter  analysis  needs  to  include 
competition  between  the  two  pathways  in  order  to  achieve  the  most  accurate  information.  Such  a 
competitive  analysis  can  be  accomplished  using  statistical  tools  well-developed  in  our 

TO 

laboratory,  but  requires  knowledge  of  the  molecular  parameters  of  the  products  and  the  barrier 
along  the  potential  energy  surface.  Again  this  points  to  have  a  more  complete  knowledge  of  the 
complete  potential  energy  surface  by  examining  this  system  computationally.  A  combined 
experimental/computational  exploration  of  this  potential  energy  surface  should  prove  insightful. 


Fig.  5 

Gd02+ 
energy 
axis)  and 
frames. 


Cross  sections  for  reaction  of 
+  Xe  as  a  function  of  kinetic 
in  the  center-of-mass  (lower  x- 
laboratory  (upper  x-axis) 
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